extern doublereal dlapy2_(doublereal *x, doublereal *y);
